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Abstract 

We present a systematic evaluation of the shape of the neutrino energy spec- 
trum produced by beta-decay of ^B. We place special emphasis on determining 
the range of uncertainties permitted by existing laboratory data and theoret- 
ical ingredients (such as forbidden and radiative corrections). We review and 
compare the available experimental data on the ^B(/3"'")^Be(2a) decay chain. 
We analyze the theoretical and experimental uncertainties quantitatively. We 
give a numerical representation of the best-fit (standard-model) neutrino spec- 
trum, as well as two extreme deviations from the standard spectrum that 
represent the total (experimental and theoretical) effective ±3(7 deviations. 
Solar neutrino experiments that are currently being developed will be able 
to measure the shape of the ^B neutrino spectrum above about 5 MeV. An 
observed distortion of the ®B solar neutrino spectrum outside the range given 
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in the present work could be considered as evidence, at an effective signif- 
icance level greater than three standard deviations, for physics beyond the 
standard electroweak model. We use the most recent available experimental 
data on the Gamow-Teller strengths in the A = 37 system to calculate the 
neutrino absorption cross section on chlorine: ac\ = (l.MibO.ll) x 10~^^ cm^ 
(±3cr errors). The chlorine cross section is also given as a function of the 
neutrino energy. The neutrino absorption cross section in gallium is 
(7Ga = (2.46if;i) X 10"^^ cm2 (±3a errors). 



PACS number(s): 96.60.Kx, 23.40.Bw, 23.60.+e, 27.20.+n 
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Typeset using REVT^ 



I. INTRODUCTION 



The solar neutrino spectrum is currently being explored by four underground experi- 
ments: the pioneering Homestake (chlorine) detector , the Kamiokande (water-Cherenkov) 
detector [|], and two gallium detectors, GALLEX g and SAGE §. 

These first-generation experiments have shown that the observed solar neutrino event 
rates are lower than expected [§, giving rise to "solar neutrino problems" ^ that cannot 
be solved within the standard experimental and theoretical understanding of the physics of 
the sun and of the electroweak interactions. 

Bahcall has shown that the ^B(/3"'")^Be (allowed) decay produces the same shape for 
the (^B) neutrino spectrum, up to gravitational redshift corrections of O{10~^), independent 
of whether the neutrinos are created in a terrestrial laboratory or in the center of the sun. 
Thus, experimental evidence for a deviation of the ^B solar neutrino spectrum from the 
laboratory shape would constitute evidence for physics beyond the standard model of the 
electroweak interactions. Indeed, powerful new experiments, such as the SuperKamiokande 
P), the SNO 0, and the ICARUS ||lOl detectors, will measure the shape of the high-energy 
{Ely ^ 5 MeV) part of the solar neutrino spectrum, which originates from the beta-decay of 
^B produced in the sun. 

The calculation of the ^B solar neutrino spectrum dates back to 1964, when it was pointed 
out [|lT| that the usual /?-decay allowed spectrum should be averaged over the intermediate 
2+ states of ®Be, as derived by the subsequent alpha-decay ^Be(2a). Experimental evidence 
for this smearing was found 23 years later in the associated positron spectrum [jl2[. The 
calculation of the spectrum has been continually improved by Bahcall and collaborators in 
ll|,13jl^, to which we refer the conscientious reader for all the details not reported here. 



The independent calculation of the B neutrino spectrum by J. Napolitano et al. [13] also 



compares well with the results in |14 



In this work, the evaluation of the *B solar neutrino spectrum is further improved, using 
recently available experimental data and new theoretical calculations. Moreover, the max- 
imum uncertainties (±3 effective standard deviations) that can affect the best estimated 
spectrum are determined. The inferred spectral shape is presented numerically and graphi- 
cally, in forms useful for fits to experimental results and for phenomenological analyses. The 
neutrino spectra, together with new data on the Gamow-Teller strengths in the A = 37 
system, are used to improve the calculations of the '^B neutrino absorption cross sections for 
chlorine and for gallium. 

The paper is organized as follows. In Sec. II, the available experimental data on the 
^B(/3"'")^Be(2a) decay chain are reviewed and discussed, in order to extract an optimal data 
set and in order to evaluate the experimental uncertainties. In Sec. Ill the best neutrino 
spectrum is calculated, including the appropriate radiative corrections. The total (theoreti- 
cal and experimental) uncertainties are used to calculate a "-|-3(t" and a "— 3(t" spectrum, 
characterizing the maximum deviations from the optimal spectrum. In Sec. IV these spectra 
are applied to a refined calculation of the *B neutrino absorption cross sections for ^^Cl and 
for ''^Ga. A summary of the work is presented in Sec. V. 
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II. EXPERIMENTAL DATA ON THE ^B{(3+fBe{2a) DECAY CHAIN 

In this section, we present a compilation of the available experimental data on beta 
and alpha decay of ®B. Then the various a-decay data are compared and used to fit the 
/3-decay data. The best-fit alpha spectrum and its ±3cr range emerge naturally from this 
overconstrained comparison. 

The /3-decay of populates 2"*" states in *Be, which then breaks up into two a-particles 



I5| . The energy levels are shown in Fig. 1. The /?, a and correlated P-a spectra have been 
investigated by several experimental groups. In addition to the solar neutrino application, 
this decay chain is of special interest because the A = 8 nuclear isotriplet (Li, Be, B) can 
be used to search for violations of the conserved vector current (CVC) hypothesis or for the 
existence of second-class currents (SCC) (see, e.g., the review |]T^), and is also an interesting 
subject for i?-matrix analyses |1T7|-[T9[|. 

For the determination of the ^B (solar) neutrino spectrum, a single precise determination 
of the alpha spectrum is, in principle, all that is required. In practice, we take advantage 
of the redundancy of the available experimental information to estimate both a preferred 
alpha spectrum and the possible uncertainties. 



A. Beta decay data 

The most recent determination of ^B /9-decay spectrum is reported (graphically) in [IT2 



The original number of detected counts in each of the 33 channels used (for momentum 
p ^ 9 MeV) have been made available by one of us (J.N.). We discuss below in detail the 
energy calibration of the f3 decay data of [0, for its particular relevance to the present 
work. 



The calibration procedure used in ||T2[ was similar to that in [^, which employed the 
same spectrograph to determine the weak magnetism correction to the ^^B /5-spectrum shape. 
In []T^ , the ^^B spectrum itself was used to set the calibration. The ^^B source was produced 
through the ^^B{d, p) reaction, with the spectrograph fixed at the same field setting used 
for the ^B data. The momentum measurement was in channels corresponding to a position 
along the focal plane. Prior studies of the spectrograph ||21| demonstrated good linearity 



between momentum and position. The ^^B raw data (~ 1.1 x 10^ counts) were collected 
in 36 momentum bins. The data were fitted with the standard allowed (3-decaj spectrum. 



along with the known recoil order corrections (see ||20| and references therein). Both the 
normalization and the offset in the position relative to the radius of curvature were left 
free in the fit. The simultaneous minimization of with respect to these two parameters 
yielded a good fit (x^/^df = 1-1), and showed no evidence for a systematic deviation in the 
residuals. The error in the absolute momentum calibration was estimated to be less than 
0.2% at la. Our conservative estimate for the maximum (±3(t) uncertainty in the energy 
calibration of the /5-decay data in [0 is 6E/^ = ±0.090 MeV. This value corresponds to 
three times a ±0.2% error at the endpoint energy, Ep ~ 15 MeV. 

As we will see, the /3-decay spectrum data of play a fundamental role in constraining 



the uncertainties of the neutrino energy spectrum. Attempts to confirm the positron spec- 
trum data, including measurements at lower momenta, would be helpful. Unfortunately, 
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existing additional data on the /3-decay as reported (only graphically) in f3-a correlation 
searches [|^,|^ are too sparse to be useful for our purposes. 



B. Alpha decay data 

Figure 2 shows our compilation of the experimental data for the delayed a-spectrum. The 
measurements with the highest statistics have been performed by Wilkinson and Alburger 



2^ , using both a thick and a thin catcher (WAl and WA2 in Fig. 2, respectively). For our 
purposes, it is sufficient to know that the typical energy loss of a-particles is ~ 100 (~ 50) 
keV in the thick (thin) catcher. The WAl data (~ 2.1 x 10® counts) have been reported by 
Warburton together with the proper channel-energy calibration formula. Barker |1^ 
has described similarly the WA2 data (~ 2.5 x 10® counts). An older measurement of the 



a spectrum was performed by Farmer and Class |25|, although reported only in a graphical 
form (Fig. 2 in |^). The spectrum labeled "FC" in Fig. 2 corresponds to a digitized form of 
their data (~ 0.5 x 10® counts). A high-statistics data set has been also collected by one of 
us (L.DeB.) and D. Wright, in the course of a recent experimental search for CVC violation 
and sec effects in the A = 8 multiplet p6[. These data (~ 1.6 x 10® counts) are labeled as 



"DBW" in Fig. 2. 

In all of the a-decay measurements we use, the experimentalists devoted considerable 
attention to the a-energy calibration, and in particular to the accurate estimate of the 
energy loss in the target. The spectra were corrected by the experimentalists for calibration 



and energy loss effects. We do not use the a-spectrum measured by Clark et al. in p7 
since data cannot be extracted from their Fig. 2 with sufficient precision to be useful for our 
purposes. 

The spectra in Fig. 2 are peaked at Ep^ak — 1.5 MeV and decrease very rapidly for 
Ea 7^ Epeak (uoticc the logarithmic scale). At LSIq!? MeV, the spectral values are decreased 
by a factor ~ 10. The interval 1.5lo!7 MeV contains ~ 90% of the experimental counts 
for each data set. Therefore the "tails" of the spectra {E^ ^ 0.8 MeV and E^ ^ 3.3 MeV 
contribute only ~ 10% to the smearing of the endpoint energies in the calculation of the 
neutrino or positron spectra from ^B decay. 



C. Discussion of experimental uncertainties 

A close inspection of Fig. 2 reveals that the four different a-decay spectra WAl, WA2, 
FC and DBW have slightly displaced peaks, with a total spread in the peak energies of about 
±0.08 MeV. Indeed, the largest uncertainty in these spectra can be ascribed to a possible 
bias h in the measured alpha-particle energy: Ea ^ E^ + b. If the intermediate 2"*^ state at 
3.04 MeV (see Fig. 1) were a narrow resonance, then b could be taken as a constant bias. In 
general, b may depend on E^, since in fact the intermediate state is not very narrow; thus 
the bias might assume slightly different values at the peak or in the tails of the a-decay 
spectra. However, in the calculation of the *B neutrino spectrum the a-spectrum tails are 
much less important than the peak region (0.8 ^ -Eq ~ 3.3 MeV), so that b is assumed to 
be a constant (6 ~ bpeak) for our purposes. 
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For the spectra WAl and WA2, Barker has given in a thorough discussion and 
an estimate of the possible contributions to b, including uncertainties due to energy loss, 
calibration, and finite resolution. A succint summary of Barker's analysis is that |6| ^ 0.05 
MeV at la. In particular, for the spectrum WA2 (thin catcher), two possible channel- 



energy calibration formulae are presented in |T9| [see Eq. (1) and (15) therein], which differ 



by ~ 0.05 MeV in the alpha-energy peak. The first calibration [Eq. (1) in Barker's paper] 
has been used in connection with the spectrum WA2 shown in Fig. 2; the second calibration 
will be used below for a further check of the sensitivity of the inferred spectrum to possible 
systematic uncertainties. For the DBW spectrum, the uncertainty in is estimated to be 
of comparable magnitude (~ 0.04 MeV) [^. The uncertainty is larger {\b\ ^ 0.1 MeV) for 
the older FC spectrum ||25|, in part as a result of the necessity of converting to graphical 
data. 

Each of the four a-particle data sets can be used to estimate the theoretical positron 
spectrum in ^B(/?"'")*Be decay. The ingredients of this calculation are the same as for the 
neutrino spectrum, except for the radiative corrections. In fact, in beta-decay the radiative 
corrections take a different form according to whether the /3-particle, or the neutrino, is 
detected. In the first case, they have been computed by Sirlin p^. In the second case. 



they have been recently evaluated by Batkin and Sundaresan [123], and will be discussed 



later (in Sec. III). The computed positron spectrum, including radiative corrections, is then 



compared to the experimental /3-decay data (33 points) as in [T^. The total number of 
counts collected (~ 0.3 x 10^) is used for normalization, reducing the number of degrees of 
freedom (Ndf) to 32. A normalized Xn (Xn = X^/^df) is then calculated for each input 
a spectrum, and this exercise is repeated also by shifting the experimental a-energy values 
{Ea Ea + b). In the calculation of Xn include the statistical errors of the /5-decay data 
but exclude, in first approximation, their energy calibration uncertainty SE/s (Sec. II. A). 
The effect of this additional uncertainty is discussed at the end of this section. 

Figure 3 shows, for each of the four measured a-particle spectra, the normalized 
to the measured positron spectrum as a function of the assumed a-particle energy bias. 
For zero energy bias (vertical dashed line) the FC, DBW, WA2, and WAl data provide 
increasingly good fits to the beta-decay spectrum. When an allowance for a non-zero bias 
b is made, all the four alpha-spectra provide equally good fits {x% — 1); modulo the shift 
Ea ^ Ea + b. The dotted curve, which also provides a good fit for zero bias, corresponds to 



the spectrum WA2 using the second Barker's calibration [0. The main difference among 



the four a-particle data sets can thus be ascribed to small biases in the measured a-energies. 

The natural choice for the optimal a-particle energy spectrum is seen to be, in Fig. 3, 
WAl with b = +0.025 MeV. This is our choice for the experimental input data in the 
calculation of the "best" neutrino spectrum. An additional variation of b equal to ±0.035 
MeV produces a Ax^ = 9 increase in the fit to the spectrum, defining a ±3cr range for b. 

We have studied the sensitivity of the x^-fit to the low-statistics bins by excluding up to 
10 bins in the high-energy part of the positron spectrum; the highest-energy bins represent 
4.5% of the experimental counts. The central values of the bias b in these fits (excluding 
some or all of the 10 highest-energy bins) are spread by ±0.005 MeV around b = 0.025 MeV, 
with a 3a error that can be as large as 0.047 MeV when all the 10 bins are excluded. 

The Kolmogorov-Smirnov (K-S) test provides a non-parametric (bin-free) way of deter- 
mining the goodness of fit of two distributions. We have therefore applied a K-S test to the 
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best-fit (unbinned) normalized spectrum. We obtain a So" error of ±0.056 MeV for the bias 
b. We adopt this error, ±0.056 MeV, as a conservative but relevant 3a estimate. 

Figure 4 shows the experimental beta-decay spectrum |]I2| (dots with la statistical error 
bars), the best-fit theoretical spectrum (solid curve, obtained by using WAl data with 
b = 0.025 MeV), and the "±3cr" theoretical spectra (dotted curves, WAl data with b = 
0.025 ±0.056 MeV). 

As discussed in Sec. II. A, the /5-decay reference spectrum is affected by a maximum 
energy calibration uncertainty 6Ej3 = ±0.090 MeV (3cr). An error ±6Ef^ corresponds to an 
error ^26Ea in the ^B(/5"'")^Be(2a) decay chain. As a consequence, the total range of the 
a-energy bias b gets slightly enlarged: b = 0.025 ±0.056 ±0.045 = 0.025 ±0.072 MeV, where 
the two (independent) errors have been added in quadrature. 

The estimated ±3a range for b (±0.072 MeV) is approximately equal to the spread in 
the peaks of the experimental alpha-decay spectra (±0.08 MeV) discussed at the beginning 
of this section. 



III. THE STANDARD NEUTRINO SPECTRUM AND ITS UNCERTAINTIES 

In this section, the main results of the present paper are presented: a best (standard) 
neutrino spectrum X{E^), together with two supplementary spectra, X'^{E^) and X~{E^), 
obtained by stretching the total uncertainties to their ±3cr limits. These spectra are given 
both in figures and in tables. The effects of the individual experimental and theoretical 
uncertainties are also discussed and illustrated graphically. 



A. Optimal neutrino spectrum and its 3(7 deviations 



The successive calculations of the ®B normalized neutrino spectrum X{E, 



-the intermediate state smearing 



besides the phase space factor 
function |Tl|,p!3|, and the forbidden corrections to the allowed transition 



included — 
the proper Fermi 



Napolitano et al. |]T2| pointed out the potential relevance of radiative corrections, al- 
though only those corresponding to the ^B positron spectrum were known at the time. 
Here we include the appropriate radiative corrections to the ^B neutrino spectrum that have 
been recently calculated in [^]. These corrections are smaller (due to a cancellation between 
real and virtual photon contributions ||2^), and with a milder energy dependence, than the 
corrections that apply when the charged lepton is detected in /3-decay [^. As we shall see 
below, their inclusion makes no significant difference in the calculation. 

The optimal input alpha-decay data, as discussed in the previous section, are taken as 
WAl with an energy bias b having a central value of 0.025 MeV and a ±3(T experimental 
uncertainty A6 = ±0.072 MeV (determined by the fit to the /3-decay data). We anticipate 
(see below) that the inclusion of the maximal theoretical uncertainties can be mimicked by 
enlarging the above range to A6 = ±0.104 MeV. The "3ct different" neutrino spectra, 
and A~, are calculated for the extreme values of A6 (±0.104 and —0.104 MeV, respectively). 

The results for the normalized neutrino spectra are shown in Fig. 5. (The spectra all 
happen to be almost coincident at about half of the maximum energy.) Numerical values of 
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A, A"*" and A are reported in Table I. A computer-readable version of Table I is available at 
the WWW site: |http:// www .sns.ias.edu/ j nb /neutrino . lit ml| . 

An alternative representation of the above neutrino spectra, in which the trivial part 
of the energy dependence is factorized out, is shown in Fig. 6 (a Kurie plot Q), where the 
ordinate is ^JXjEy. Notice the deviation from a straight line; the deviation is primarily due 
to the smearing over the intermediate broad state of ^Be. 

A representation of the integral spectrum is shown in Fig. 7, where the fraction / of 
'^B neutrinos produced above a given neutrino energy threshold Et^ \j = lE^t^dE^X^E^)] is 
plotted as function of Ef^. 



B. Spectrum uncertainties: Experimental and theoretical components 

The effect of varying the input alpha-decay data with respect to the optimal choice (WAl 
with bias b = 0.025) is shown in Fig. 8, where the solid line represents the standard spectrum 
A, and crosses are placed at representative points along the spectra A obtained with WAl, 
WA2, FC and DBW input data (with no bias). The differences AA = A — A are very small, 
and can be best appreciated in an expanded scale in Fig. 9, where the dimensionless quantity 
AA/Apeafc is plotted [Xpeak = maxA(E^)]. 

Figure 9 also shows the ±3a deviations (A^ — X)/ Xpeak- These deviations are similar 
to sinusoidal curves, with a maximum amplitude of ~ 2.5%. The average value of the 
absolute deviation is thus (|AA|/Apeafc) — (2/7r)|AA|maa;/Apeafc — 1-6% at 3a. The difference 
between the Bahcall and Holstein fl^ spectrum Abh and the best-fit spectrum A can also 
be represented well by a sinusoidal curve (like those shown in Fig. 9). The amplitude of the 
difference Abh — A is ~ 0.7a, to be compared with the effective 3a differences A^ — A shown 
in Fig. 9. Similarly, the difference between the spectrum An calculated by J. Napolitano et 
al. in [0 and the standard spectrum A is about l.Aa. 



Notice that the curve FC in Fig. 9 is somewhat irregular and "out of phase" (about | 
of a semiperiod) with respect to the ±3cr sinusoids. The irregularity is due to the scatter of 
the FC data points. The dephasing can be traced back to the fact that the FC data set has 
very few low-energy counts, being limited to E^ ^1.2 MeV. We have verified that the curve 
FC in Fig. 9 is "rephased" if the FC alpha spectrum of Fig. 2 is artificially prolonged to 
lower energies. We have also computed neutrino spectra with "mixed alpha data," namely 
with WAl data for E^ > Ep^ak matched to DBW data for E^ < Ep^ak (and viceversa). 
We found that the spectral differences can always be reabsorbed in a uniform bias h to 
a very good approximation, with residual differences of ^ 0.015 MeV with respect to the 
"unmixed" original data. We conclude that the low-energy part of the experimental alpha- 
decay spectrum (that affects particularly the high-energy tail of the neutrino spectrum) is 
sufficiently well known for our purposes, and that the alpha spectra uncertainties can be 
effectively parametrized as a uniform energy offset h. 



We have not divided the spectrum by the Fermi function as is usual in plotting beta-decay 
spectra. In the present case, the Fermi function would have to be averaged over a range of positron 
energies because of the width of the final (2+) state in ^Be. 
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The effects of radical assumptions about the correctness of the theoretical calculations 
beyond the (phase space) x (Fermi function) approximation are shown in Fig. 10, on the same 
scale as Fig. 9. The curves labeled 1 and 2 represent the shifts AA/Apeafc obtained by setting 
to zero the forbidden or the radiative corrections, respectively. For curve 3, the radiative 
corrections were (inappropriately) assumed to be the same as for the positron detection 
case (as was also done in [0); the present exercise is intended to account roughly for a 
hypothetical situation in which the cancellation between real and virtual photon contribution 
might not be as effective as computed in ||29|| . 

The maximum theoretical deviation (curve 1 in Fig. 10) is obtained by excluding the 
forbidden corrections altogether. We have verified that this deviation can be mimicked by 
recalculating the neutrino energy spectrum with an additional energy bias, A6theo — 0.075 
MeV. This value is comparable to the estimated 3a experimental uncertainty evaluated 
in Sec. II C, A6exp — 0.072. We assume that the maximum theoretical offset, Afetheo, 
corresponds to an "effective" 3a statistical significance. Our final best-estimate for the bias 
b to be applied to the reference WAl a-decay spectrum is therefore: b = 0.025 MeV ± 
(A62^p + A6t\gj5 = 0.025 ± 0.104 MeV, as anticipated in Sec. III.A. 

The assignment of a "Sex level of significance" to the offset that parameterizes the total 
effect of the forbidden corrections is a plausible estimate. However, there is no rigorous way 
to estimate a 3a theoretical uncertainty. Our estimate is motivated by the fact that the 
calculation of forbidden corrections in beta-decay is not made purely from first principles 
(strong and electroweak lagrangian), unlike the QED radiative corrections. The evaluation 
of the forbidden terms makes use of approximate symmetries to expand the weak nuclear 
current in terms of a few nuclear form factors, that are evaluated using nuclear models and 
difficult f3-a correlation experiments (see |14] and references therein). 

If the reader prefers to adopt a smaller or larger estimate for the theoretical uncertainty, 
a simple prescription for estimating the changes in the inferred uncertainties is given below. 
The prescription is based upon linear error propagation and generally gives agreement with 
an exact numerical calculation to a fractional accuracy of about 5% or better in the change 
induced by the rescaling of the error. Recalculate the total 3a range for the bias, Ab = 
(A6g^p + A6^[^gQ)2, with A6exp = 0.072 MeV and the reader's preferred estimate for Afotheo- 
Then rescale by a factor / = A6/(0.104 MeV) all the ^B-related 3a total uncertainties 
quoted in this paper. For example, the 3a spectral deviations from the best-estimated 
neutrino spectrum shape, A^(-E,^) — \{E^), become /[A^(£'i,) — X{E^)]. Notice that the 
rescaling factor is at least / = 0.69, which is obtained by setting Afetheo = 0. 

A final remark is in order. The relative contribution of the different 2+ states of ^Be 
in the a-decay spectrum (see Fig. 1) has been analyzed by Barker [|T3,|r^ and Warburton 
||T8| within the i?-matrix formalism. Their results are not in complete agreement, the fitted 
amplitude of the intermediate states being sensitive to the input data (see, e.g., the discussion 
in [0). In particular, i?-matrix fits are very sensitive to the absolute energy calibration, 
as well as to the tails of the alpha-decay spectrum. In our calculation of the neutrino 
spectrum, the absolute a-energy is allowed to vary within the quoted uncertainties (±0.104 
MeV). The tails of the alpha decay spectrum are not decisive for our purposes. We are not 
really interested in separating the relative contributions of the 2+ states; we just rely upon 
the global population of the 2"*" states as derived from the experimental a-decay data. We 
conclude that the disagreement between the theoretical i?-matrix fits p!7HT9|] is not an issue 
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in our calculation of the ®B neutrino spectrum. 



IV. ABSORPTION CROSS SECTION OF NEUTRINOS IN ^^Cl AND ^^Ga 



In this section, we present improved calculations of the *B neutrino absorption cross 
section for chlorine and for gallium. Recent calculations of the chlorine absorption cross 
section were made by Bahcall and Holstein |T^, Garcia et al. and Aufderheide et 

al. |31[]; the results of earlier calculations can be found in [lTT| , p!3| , |32| -[36|l . The most recent 



previous calculation of the gallium cross section crca was by Bahcall and Ulrich . 

Transitions to excited states dominate the total cross section in either of the absorp- 
tion processes '^^Cl(z/, e)^^Ar and ^^Ga(z/, e)^^Ge. The Gamow-Teller transition strengths, 
i?(GT), can be estimated from the rates of the analogous charge-exchange (p, n) reactions. 
For the A = 37 system, these transition matrix elements can be determined experimen- 
tally by studying the ^^Ca(/5"'")'^^K transition [^,|TT|, which is the isospin mirror process 

^ for a more 



of ^^Cl(z/, e)^^Ar. The interested reader is referred to the recent review in 
extensive discussion of these processes. 



A. Absorption Cross Section for Chlorine 



Including for the first time forbidden corrections, Bahcall and Holstein calculated 
the ^B neutrino cross section on chlorine and obtained: 



aci = (1.06 ± 0.10) X 10"^^ cm 



(1) 



The quoted uncertainties represented the maximum estimated error (Scr). The calculation 
made use of the -B(GT) values derived from the ^''Ca /3-decay, as reported by Sextro et al. 
in [Q. The estimated Sex error (±0.10) had two components, ±0.08 from ^Be a-decay data 
and ±0.06 from ^''Ca /3-decay data uncertainties, to be added in quadrature. 

T^ , and the spectra A, A"*", and A~ reported in 
The best-estimate calculated cross sections 



Using the same low-energy data |Q as in 
Table I, we find aci = (1.08 ±0.15) x lO'^^ 
differ by 2%. The 3o" error component from ®Be(2a) decay data is ±0.07 



Recently, new experiments have been carried out to determine more precisely the i?(GT) 



strengths in ^^Cl(p, n)^^Ar and '^^Ca(/3"'")^'^K 0] processes. Taken at face values, the 
i?(GT) strengths derived by the different experiments [^,0,^ were not in good agreement. 
Critical examinations ||39| of the data analyses, as well as supplementary data , may have 
led to a satisfactory understanding of the low-energy levels and their i?(GT) strengths 
in the A = 37 system. 



Using the latest available data |3l|] and the neutrino spectra A*^^-* calculated here, we find 



(Tci = (1.14 ±0.11) X 10 



-42 



cm 



(2) 



Equation (2) represents our best estimate, and the associated ?>a uncertainties, for the ^B 
neutrino absorption cross section on chlorine. The contribution to the total error from the 
measured i?(GT) values is assumed to be ±0.08, as in the analysis The difference 
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between the values of the chlorine absorption cross section in Eq. (1) and Eq. (2) is 7%, well 
within the quoted errors. 

Figure 11 shows the energy dependence of our our best-estimate chlorine cross section 
(solid line). Values of cxci for representative neutrino energies are also given in Table II. The 
difference between the present and the previous calculation of the energy- dependent cross 
section by Bahcall and Ulrich |^ (dashed line in Fig. 11) is less than 20% for Ei, < 16 MeV. 
The differences are largest at the highest energies since the newer data include transitions to 
higher excitation states in ^^Ar that were not determined in the previous experiments (see 
|3g,^g). 



To give the reader some perspective on how the B neutrino spectrum and the chlorine 
absorption cross section have changed with time, we give in Table III all the published 
values of cxci with which we are familiar. The calculated cross sections have been approxi- 
mately constant, within the estimated errors, since 1978, although there have been numerous 
refinements (which are described in |P!B,H, SB| , |SD| , PT| ). The reasons for the relatively signif- 



icant change in the 1978 best-estimated value [jT3[ with respect to the earlier calculations 
Tl| , |33| -p5| are described in the last paragraph of Sec. IV. B. 3 in [|1^ . 



B. Absorption Cross Section for Gallium 

The ®B neutrino absorption cross section for gallium that is widely used was calculated 
by Bahcall and Ulrich and is: 

cTGa = (2.431?:}) X 10-^2 cm^ , (3) 

where the quoted uncertainties represented the maximum estimated errors (3(t). The B{GT) 
values used in the quoted calculation were taken from a ''^Ga(p, n)^^Ge experiment per- 
formed by Krofcheck et al. ^T|. 

The only important recently-published experimental development with which we are 



familiar is the recent ^^Cr source experiment for the GALLEX detector [42|. Rata and 
Haxton have shown that the measurements with the chromium source are consistent 
with the i?(GT) values for the first two excited states that were inferred by Krofcheck et al. 
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Therefore we repeat the Bahcall-Ulrich calculation [Eq. (3)] using the best ^B neutrino 
spectrum from the present paper. We find: 

(TGa = (2.46^?:}) X 10-^2 ^^2 (4) 

The change in the best-estimate cross section is only ~ 1% [relative to Eq. (3)], which is much 
smaller than the guessed systematic errors, that represent uncertainties in the interpretation 
of the (p, n) measurements. 



V. SUMMARY 

In the previous sections, the spectrum of neutrinos produced in the ^B(/?"'')*Be(2a) decay 
has been computed, using state-of-the-art theory of beta-decay. The laboratory data on the 
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associated positron spectrum have been used to choose an optimal data set among the 
different measured ^Be(2a) decay spectra. The experimental and theoretical uncertainties 
can both be represented well as an energy shift (b) in the a-decay data. The total ±3a 
range for this shift (bias) has been conservatively estimated to be ±0.104 MeV. A best- 
fitting standard spectrum, A, has been computed, as well as the "effective ±3cr" neutrino 
spectral shapes A"^ and A^ [A"^ — A = 3 standard deviations]. The standard spectrum A 
differs by about 0.7cr from the Bahcall and Holstein neutrino spectrum []T1| and by about 
1.4(7 from the Napolitano et al. spectrum [0. 

The neutrino absorption cross section for chlorine calculated with the best-fitting 
spectrum derived here and with the most recent data on the low-lying states in the A = 37 
system is aci = (1.14 ± 0.11) x 10"^^ (3a). This result is in agreement with the 
estimates derived in 1964 (see Table III). The best-estimate gallium absorption cross section 
is aca = (2.46;?;}) x lO'^^ ^^^2 ^^^y 

Many readers would prefer to quote la rather than 3(j errors. To a good approximation, 
la errors can be obtained from our quoted 3cr values by dividing by three. Moreover, 
[A±(E,) - A(E,)]i. ~ [A±(E,) - A(E,)]3./3. 

All the available experimental data on the ^B(/5+)^Be(2a) decay are consistent within 
the quoted uncertainties. Higher order contributions in the theoretical calculation of the ^B 
neutrino spectrum should be very small. Any measured deviation of the ^B solar neutrino 
spectrum in excess of the conservative limits given in this paper could be considered as 
evidence for new physics beyond the standard electroweak model, at an effective significance 
level greater than three standard deviations. 
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TABLES 



TABLE L The spectrum A of solar neutrinos from the decay of ^B, together with the spectra 
associated with the total itScj uncertainties. The neutrino energy is MeV and \'^^\Ey) is the 
probability that a neutrino with energy Ey is emitted between £'^^±0.05 MeV. A computer-readable 
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TABLE II. Values of the absorption cross section in chlorine, in units of 10"^^ cm^, for repre- 
sentative values of the neutrino energy. The second column refers to the calculation in the present 
paper. The third column refers to the Bahcall-Ulrich (BU) ||3^ calculation. 



Ey (MeV) 




o-ci(BU) 


1 


5.21E+00 


5.21E+00 


2 


3.70E+01 


3.70E-I-01 


3 


1.02E+02 


1.15E-I-02 


4 


2.23E+02 


2.63E-I-02 


5 


5.38E+02 


5.63E+02 


6 


1.44E+03 


1.52E+03 


7 


4.62E+03 


4.76E+03 


8 


l.OlE+04 


1.02E-I-04 


9 


1.85E+04 


1.79E+04 


10 


3.00E+04 


2.77E-I-04 


11 


4.45E+04 


3.97E+04 


12 


6.21E+04 


5.38E+04 


13 


8.27E+04 


7.00E+04 


14 


1.06E+05 


8.83E+04 


15 


1.33E+05 


1.09E+05 


16 


1.62E+05 


1.31E+05 


18 


2.28E+05 


1.81E+05 


20 


3.05E+05 


2.38E+05 


30 


8.20E+05 


6.11E+05 



TABLE III. Values of the ^B neutrino absorption cross section for chlorine (o"ci), as calculated 
by various authors. The first and second (when given) error components, ecT and ee, are to 
be added in quadrature; they refer to the estimated uncertainties from the Gamow-Teller (GT) 
strengths and from the ^B neutrino spectrum, respectively. When the definitions of the errors 
given in the original papers were sufficiently precise, we have indicated (in parentheses) that we 
are quoting 3(7 errors. 



Year 



Author(s) 



Ref. 



fjci ± ecT ± eB (10 cm^) 



1964 
1964 
1966 
1974 
1977 
1978 
1981 
1986 
1991 
1994 
1996 



Bahcall 

Bahcall 

Bahcall 

Sextro et al. 

Haxton and Donnelly 

Bahcall 

Rapaport et al. 
Bahcall and Holstein 
Garcia et al. 
Aufderheide et al. 
Bahcall et al. 







m 
m 



m 



m 




m 



m 

This work 



1.27 ±0.31 
1.30 ±0.29 
1.35 ±0.10 
1.31 

1.27 ±0.22 ±0.06 

1.08 ±0.10 
0.98 ±0.07 

1.06 ± 0.06 ± 0.08 [3cr] 

1.09 ±0.03 
1.11 ±0.08 [3o-] 

1.14 ±0.08 ±0.08 [3o-] 
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FIGURES 



FIG. 1. Energy levels in the ^B(/3^)®Be(2a) decay chain (not to scale). 

FIG. 2. Compilation of ®Be(2a) decay data. The bin widths are different for different experi- 
ments. The data WAl and WA2 are shifted on the vertical axis. 

FIG. 3. Values of the normalized chi-square in a fit to the experimental positron spectrum, 
using the input alpha-decay data of Fig. 2, with an allowance for a possible bias, b, in the detected 
a-particle energy. The curves are remarkably similar, modulo a constant bias. 

FIG. 4. Experimental data on the positron spectrum, together with the best fit and the ±3(7 
fits, corresponding to WAl alpha-decay data within the bias range b = 0.025 ± 0.056 MeV. 

FIG. 5. The best-estimate (standard) neutrino spectrum A, together with the spectra 
allowed by the maximum (±3(t) theoretical and experimental uncertainties. 

FIG. 6. The spectra A, A^, and A^, shown as a Kurie plot. 

FIG. 7. Fraction of ^B neutrinos produced with energy Ei, above a given threshold Eti^^. 

FIG. 8. Variations in the standard neutrino spectrum A (solid line) induced by different input 
data sets (crosses). 

FIG. 9. Variations A A in the standard neutrino spectrum A induced by different input data 
sets, divided by the peak value of A (Xpeak)- The maximum (±3cr) differences A^ — A are also 
shown. 

FIG. 10. Variations AA in the neutrino spectrumy induced by drastic changes in the theoretical 
computation, shown in the same scale as Fig. 9. See the text for details. 

FIG. 11. Absorption cross section in chlorine (solid line) as a function of the neutrino energy. 
The dashed line refers to the Bahcall-Ulrich calculation. 
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